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Abstract —In this paper, we present a decision level fused 
local Morphological Pattern Spectrum(PS) and Local Binary 
Pattern (LBP) approach for an efficient shape representation 
and classification. This method makes use of Earth Movers 
Distance(EMD) as the measure in feature matching and shape 
retrieval process. The proposed approach has three major phases 
: Feature Extraction, Construction of hybrid spectrum knowledge 
base and Classification. In the first phase, feature extraction of the 
shape is done using pattern spectrum and local binary pattern 
method. In the second phase, the histograms of both pattern 
spectrum and local binary pattern are fused and stored in the 
knowledge base. In the third phase, the comparison and matching 
of the features, which are represented in the form of histograms, 
is done using Earth Movers Distance(EMD) as metric. The top-n 
shapes are retrieved for each query shape. The accuracy is tested 
by means of standard Bulls eye score method. The experiments 
are conducted on publicly available shape datasets like Kimia-99, 
Kimia-216 and MPEG-7. The comparative study is also provided 
with the well known approaches to exhibit the retrieval accuracy 
of the proposed approach. 

Keywords—Pattern spectra^ Earth Movers Distance^ Histogram 
matchingy Shape retrieval. Local Binary Pattern. 

1. Introduction 

In the field of computer vision, the process of under¬ 
standing the scene involves object extraction, representation or 
feature extraction, matching and classification. In particular, 
object representation and classification constitute a deeply 
entrenched and ubiquitous component of any automated in¬ 
telligent systems. In order to classify the object, representing 
the object in suitable form with dominant features is very 
much essential. It is observed from the literature that the 
shape is one such dominant feature for object representation. 
A good shape representation should be compact and should 
retain the essential characteristics of the shape. It should be 
invariant to rotation and scale. In this paper, the representation 
process takes shape of the object as feature since shape is an 
important,effective and rich feature which is used to describe 
both image boundary as well as content and hence plays a 
significant role in classification task. 

Shape representation and description techniques can be 
generally classified into two broad categories: Contour-based 
methods and Region-based methods. Contour based shape 
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techniques exploit shape boundary information, where as, in 
region based techniques, all the pixels within a shape region 
are taken into account to obtain the shape representation. 
Contour-based approaches are more popular than region-based 
approaches in literature since human beings are thought to 
discriminate shapes mainly by their contour features. However, 
region-based methods are more robust as they use the entire 
shape information. These methods can cope well with shape 
defection which arises due to missing shape part or occlusion. 
Skeleton or medial axis transformation is an example for region 
based technique. In our work, both contour as well as skeletons 
are used to take advantage of both the methods. We present 
a novel approach devised based on the morphological pattern 
spectrum to obtain pattern spectrum for every binary shape 
making use of their contour and skeleton. Then the local binary 
pattern(LBP) of binary shape is computed. The Earth Movers 
Distance(EMD) measure is used to compare and match the 
pattern spectrum of the training set and the query shape. 

The remaining part of the paper is organized as follows. In 
section 2, the review of the related works are brought down. 
In section 3, the technique of morphological patter spectrum 
is given. In section 4, an insight into local binary pattern is 
presented, followed by the description about the EMD metric 
in section 5. The proposed approach is given in section 6. In 
section 7, the experimental set-up along with the discussion of 
the results are brought out and conclusion is given in section 
8 . 

II. Review of Related work 

The general principle adopted for shape representation 
and classification involves two basic steps: representation of 
an object retaining the dominant features and matching the 
features with a suitable distance measure. We have seen several 
methods in the literature, which have been developed for shape 
representation and classification. More in particular from the 
field of mathematical morphology introduced by Matheron 
m and Serra | [20l . It is found that the medial axis transfor- 
mation(MAT) or skeleton is one such popular descriptor based 
on morphology for shape represenation and classification. The 
skeleton of an object is the locus of the centres of the maximal 
disks that can be inscribed inside the object boundary. Skeleton 


is a very useful shape descriptor and is essential for shape 
representation and analysis in many application areas such as 
content-based image retrieval systems, character recognition 
systems, circuit board inspection systems and analysis of 
biomedical images Skeletons can capture articulation of 
shapes in better form than contours |T^ . 

Skeletons can be organized in the form of attribute-relation 
graphs (ARG) for the matching purpose. One of the examples 
for ARG is shock graphs introduced by Siddiqi et al. p4| pS] . 
Shock Graph is an abstraction of skeleton of a shape onto 
a Directed Acyclic Graph (DAG). The skeleton points are 
labeled according to the radius function at each point. Shock 
graphs are constructed using the specialized grammar called 
Shock Grammar | [23| . In the skeleton, the branch points, 
end-points, and skeleton segments contain both geometrical 
and topological information. These primitives are referred as 
shocks (T7| p9| . The concept of bone graph is an extension 
of shock graph. Bone graphs retain the non-ligature structures 
of the shock graph and are more stable 0 - There are several 
algorithms which are proposed in order to have an efficient 
match of these skeletal graphs. Sebastian et al fT^ pro¬ 
posed a matching approach for node-attributed trees which 
measures the edit distance between two graphs. Since this 
method involves complex edit operation, it suffers from high 
computational complexity. 

The techniques which involve graph matching by finding 
the correspondence between nodes through the conversion of 
skeleton graphs to skeleton trees require heuristic rules to 
select the root node 0 GD- The major drawback of this 
method is that a small change in the shape causes the root 
to change that results in significant change in the topology of 
the tree representation. Apart from this, the conversion from 
graph to a tree structure results in loss of significant structural 
information and hence leads to wrong match Q. In 2008, Bai 
and Latecki jlj proposed a method based on the path similarity 
between the end points of the skeleton graphs. In this method, 
the geodesic paths between skeleton end points are obtained 
and are matched using Optimal Subsequence Bijection (OSB) 
method. Unlike other methods, this method does not convert 
the skeleton into a graph and matching the graphs,which is 
still an open problem. This approach addresses the problem of 
shapes having similar skeleton graphs, but different topological 
structures and having different skeletons with visually similar 
shapes. This approach works well even in the presence of 
articulation and contour deformation. 

Demicri et al. 0 proposed a technique called skeletal 
shape abstraction from examples. A many to many correspon¬ 
dence between the object parts is used to match the object 
skeletons. This method is invariant to shape articulation and 
appears to be more suitable for object detection in edge images. 
However, the tree abstractions are not specific enough to pre¬ 
vent hallucinating of target object in clutter. The dissimilarity 
value between the object trees is computed by performing 
many to many matches between the vertices of the trees. Earth 
Movers Distance (EMD), under LI norm, is used as a metric 
for the matching process. 

Shu and Wu p2| proposed contour points distribution 
histogram (CPDH) as the shape descriptor. Though it is not the 
best one compared to the other state-of-art approaches in terms 
of retrieval rate, it is relatively simple and has low time com¬ 


plexity. Wang et al. p7| proposed a shape descriptor which 
represents the object contour by a fixed number of sample 
point. Each sample point is associated with a height function. 
This method is capable of handling nonlinear deformation of 
the objects. Shen et al. pA] proposed a method which finds the 
common structure in a cluster of object skeleton graph. This 
method outperforms other methods and also suitable for large 
datasets. However, the time complexity of this method is high 
due to agglomerative hierarchical clustering. 

HI. Morphological Pattern Spectrum 
Representation 

The original idea of pattern spectrum proposed by 
Maragoas is based on Serra’s Mathematical morphology 
filters, that describe the distribution of local figure thickness. 
The pattern spectrum is a morphological tool that gives the 
quantitative information about the shape and sizes of the ob¬ 
jects in the image. It is invariant to the rotation and shift of the 
object and can be directly represented using histograms. The 
size distribution can be represented in the form of histogram 
for further processing. The additional attributes such as length, 
orientation and area can be included along with the shape and 
size. Pattern spectrum focuses on spatial distribution of the 
skeleton points with respect to the object contour. In order to 
construct the pattern spectrum, the shape need to be subjected 
to skeletonization. Then for each skeleton point x, the radius 
of the maximum circle which can be inscribed inside the shape 
contour with x as the centre is determined. The radius function, 
associated with each skeleton point is represented in the form 
of histogram forming the feature vectors. To illustrate, the 
skeleton of a shape obtained by applying distance transform 
method 0 and its corresponding thickness map with respect 
to the object boundary is shown in figure 0b) and (c). The 
thickness map is later transformed into 10 bin histogram as 
shown in figure 0d) 



Fig. 1. (a)Input shape, (b)Skeleton of shape (a), (c) Maximal disc inscribed 
on the skeleton of the shape(a), (d) Histogram representing the radius of discs 
and the number of such discs 


IV. Local Binary Pattern(LBP) Method 

The Local Binary Pattern (EBP) is an excellent metric 
which focuses on the spatial structure of the image GD- It is 
a powerful tool for rotation invariant texture classification. The 
EBP computation on a gray scale image is done as follows. 
Let gc be the pixel for which the EBP value need to be 
calculated. Let gp, {p = 0, ...,P — 1) correspond to the 
gray values of P equally spaced pixels on a circle of radius 
R{R> 0) that form a circularly symmetric neighbour set of 
gc . Without losing information, the gray value of the center 
pixel gc can be calculated as, 

T= {gc,go-gc,gi-gc,---,gp-9c)- 


( 1 ) 






Assuming that the differences — Qc are independent of 
gc, T can be factorized as 


measures the least amount of work needed to fill the holes 
with earth. 


T « t{gc)t{go - gc, gi- gc,---, gp -1 - gc) (2) 


The factorized distribution is only an approximation of the 
joint distribution as exact independence is not there in practice. 
Hence, though there is a small loss in information existence, 
this can be tolerable to achieve invariance with respect to shifts 
in gray scale. The invariance is achieved with respect to the 
scaling of the gray scale by considering just the signs of the 
differences instead of their exact values. 

T fa t{s{go - gc), s{gi - gc),..., - gc)) (3) 


where 


s{x) 


1, if x>0 \ 
0, if X < 0 j 


By assigning a binomial factor 2^ for each sign s{gp — gc), 
a unique transform LBPp^r is obtained that characterizes the 
spatial structure of the local image texture: 

p-i 

LBPp,n = J2s{gp-gc)2P (4) 

p=0 

LBPp^R operator is invariant against any monotonic transfor¬ 
mation of the gray scale. LBPp^p value remains constant as 
long as the order of the gray values in the image stays the 
same. 

The LBPp^p operator produces 2P different output values, 
corresponding to the 2^ different binary patterns that can be 
formed by the P pixels in the neighbour set. When the image is 
rotated, the gray values gp will correspondingly move along the 
perimeter of the circle around go. Hence, rotating a particular 
binary pattern naturally results in a different LBPp^p value. 
To make it rotation invariant, define the local binary pattern as 


To illustrate, let P = (pi,u;pi),..., (pm^'^Pm) be the first 
distribution with m clusters, where Pi is the cluster rep¬ 
resentative and wpi is the weight of the cluster and Q = 
(gi, regi),..., , re) be the second distribution with n 

clusters. In other words, p^and g^ typically represent bins in 
a fixed partitioning of the relevant region of the underlying 
feature space, and the associated wpi and wqi are a measure 
of the mass of the distribution that falls into the corresponding 
bin. Let D= dij be the ground distance matrix where dij is 
the ground distance between clusters pi and g^. We need to 
estimate the flow F = [fij], with fij be the fiow between pi 
and Qj , that minimizes the overall cost. The Minimal cost 
(WORK), between two distribution is calculated as follows: 


WORK{P, Q, F) = 

i=i j=i 

Subject to the following constraints 


E fa - 1 <i<m ; E ^ J ^ 


m 


j=l i=l 

m n / m n 

EE/*i = E^p-E^« 

i=ij=i \i=i j=i 


( 6 ) 


Constraint(l) specifies that the supplies can be moved from 
P to Q and not vice versa. Constraint(2) limits the amount 
of supplies that can be sent by the clusters in P to their 
weights. Constraint(3) limits the clusters in Q to receive no 
more supplies than their weights and constraint(4) forces to 
move the maximum amount of supplies possible. The amount 
of fiow is called as the total fiow. Once the transportation space 
problem is solved, and the optimal fiow F is found, the earth 
movers distance is defined as the resulting work normalized 
by the total fiow: 


LBPp^p = min { {LBPp^p, i) \i = 0,1,... P - 1} ; (5) 

where ROR(x, i) performs a circular bit-wise right shift 
on the P-bit number by i times. In terms of image pixels, (5) 
simply corresponds to rotating the neighbour set clockwise so 
many times that a maximal number of the most significant bits, 
starting from g^-i, is 0. 

In our work, the LBP is computed for the object pixels 
with 8 neighbours. These LBP values are represented in the 
form of histogram. 

V. Earth Movers Distance 


s;p\m \;p\n 7 n 

EMD{P, Q) = f (7) 

l^i=i 2^j=i Pj 

VI. Proposed approach 

In the proposed work, a combined approach of morpho¬ 
logical pattern spectrum and local binary pattern are used to 
extract the shape features. The LBP of the shape is computed. 
Both pattern spectrum and LBP features are represented in 
the form of histogram forming the hybrid pattern spectrum 
knowledge base. The EMD distance measure with LI ground 
distance is used for comparing the features represented in the 
form of histogram. 


The Earth Movers Distance (EMD) is defined as a minimal 
cost that must be paid to transform one distribution into the 
other It is more robust than histogram matching tech¬ 

niques and can operate on variable-length representations of 
the distribution. Eor instance, if Ha and Hb are two histograms, 
EMD is the minimum amount of work needed to transform 
histogram Ha towards Hb. Given two distributions, it can be 
seen as a mass of earth properly spread in space and the other 
as a collection of holes in the same space. Then, the EMD 


A. Feature extraction by Morphological Pattern spectrum 

The morphological pattern spectrum is computed for the 
skeleton of the shape as well as the contour. The steps involved 
in this process are given below 

1) The skeleton of the shape is extracted using skeleton 
strength map(SSM) Q which is based on the Eu¬ 
clidean distance transform method. 



2) The pattern spectrum need to be computed for each 
skeleton point. The radius function of each skeleton 
point, with respect to the shape boundary is obtained. 
The radius function value for a pixel x is the radius of 
the maximal disc inscribed inside the object with x as 
a skeleton point. The radius value for all the skeleton 
points forms the pattern spectrum of the shape. 

3) To obtain the radius values for the skeleton points, 
skeleton strength map(SSM) 0 based on the Eu¬ 
clidean distance transform is used. The forward mask 
shown in Figure [^a) and 90 degree rotated backward 
mask given in Figure [^b) are used as structuring 
elements. The forward mask is applied linearly over 
the block, which is of same size as the shape, in the 
forward direction i.e., from left to right and top to bot¬ 
tom, resulting in forward pass. The backward mask 
is applied over the image in reverse direction from 
right to left and bottom to top resulting in backward 
pass. The block is aligned with the shape skeleton 
and the values in the distance map, corresponding to 
the skeleton point, represents the radius. The result 
of forward pass and reverse pass of the block and 
the equivalent distance maps are depicted in figure 

I3d)-(f). 


The radius function is calculated using a 3 x 3 
structuring element. Results of forward and backward 
mask is shown in figure 

8) The radius function values obtained for each contour 
pixel is represented in the form of histogram forms 
the pattern spectrum knowledge base of contour 
points. 



Fig. 4. a) Block B (b) Result of the forward pass on B (c) Result of the 
backward pass on B , (d), (e), (f) represents the corresponding distance maps 
of (a), (b) and (c) 
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Fig. 2. (a) Forward mask (b) Backward Mask 
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Fig. 3. Shape B (b) Result of the forward pass on B (c) Result of the 
backward pass on B , (d), (e), (f) represents the corresponding distance maps 
of (a), (b) and (c) 




B. Feature extraction by LBP method 

The EBP value is computed for each shape pixel as follows. 
For each pixel, identify the 8 neighbours of the pixels. Obtain 
the EBP sequence of 8 neighbours. Compute the decimal 
equivalence of the this binary string. The string is subjected 
to circular right shift by 1, and then its decimal value is 
calculated. This is repeated for 8 times obtaining 8 decimal 
values. The binary string corresponding to smallest decimal 
value is taken forming the EBP value of the pixel under 
processing. The array of LBP values for each shape pixels 
is now represented in the form of histogram. This process 
is repeated for each shape in the training set forming the 
LBP knowledge base of the training set. Thus every shape 
is described by three histograms obtained due to Skeleton, 
Contour spectrum and LBP and thus forms a fused feature 
for every shape. 


C. Classification 

The Earth movers distance (EMD) is used as the distance 
measure for comparing the shape features which are repre¬ 
sented as histogram. Given a shape feature for its classification, 
the match is obtained by EMD distance as discussed in section 
V. 


Given the query shape, the classification is done as follows. 


4) The pattern spectrum thus obtained is converted to a 
histogram with 10 bins represents shape feature. 

5) The above process is repeated for all the shapes in the 
training set to form the Pattern spectrum knowledge 
base of skeleton points. 

6) A bounded box is fitted to the shape to eliminate the 
non shape portion of the image and the contour of 
the shape is extracted 

7) For each contour pixel, the pattern spectrum is com¬ 
puted with the bounding box as the reference line. 


1) Extract the skeleton pattern spectra (SPS) of the test 
shape and obtain its histogram, say St 

2) Extract the contour pattern spectra(CPS) of the test 
shape and obtain its histogram, say Ct 

3) Compute the LBP of the test shape (LBS) and obtain 
its histogram, say Lt. 

4) Compute the distance between the test shape with the 
sample shapes of the training set. Let Ds, Dc and 
Di, be the distances obtained due to EMD distance 


measure. 


















TABLE L 


Top 10 CLOSEST matching shapes 


TABLE IT 


Top 11 CLOSEST matching shapes 


Approach 

1st 

2nd 

3rd 

4th 

5th 

6th 

7th 

8th 

9th 

lOtl 

tota 

sc M 

97 

91 

88 

85 

84 

77 

75 

66 

56 

37 

756 

CPDH+EMD 
(Eucl) [g 

96 

94 

94 

87 

88 

82 

80 

70 

62 

55 

808 

CFUH+EMD 
(shift) g 

98 

94 

95 

92 

90 

88 

85 

84 

71 

52 

849 

Gen. model 

|M 

99 

97 

99 

98 

96 

96 

94 

83 

75 

48 

885 

Proposed 

Approach 

99 

97 

97 

88 

88 

86 

86 

90 

80 

77 

888 


Approach 

1st 

2nd 

3rd 

4th 

5th 

6th 

7th 

8th 

9th 

lOtl 

1 lltl 

1 Total 

sc 

214 

209 

205 

197 

191 

178 

161 

144 

131 

101 

78 

1809 

CPDH+ 
EMD 
(Eucl) ig 

214 

215 

209 

204 

200 

193 

187 

180 

168 

146 

114 

2030 

CEUH+ 

EMD 
(shift) (22] 

215 

215 

213 

205 

203 

204 

190 

180 

168 

154 

123 

2070 

Propose 

Approach 

216 

209 

205 

195 

195 

197 

188 

180 

179 

163 

152 

2079 


5) The distance matrices DS, DC and DL are fused to 
form the resultant matrix DR as follows 

DR = aDS + PDC + -fDL ( 8 ) 

where the values of a, /3 and 7 are deducted experi¬ 
mentally, which are specific to the dataset. 

6 ) The HMD distance values are arranged in the ascend¬ 
ing order of the distance to obtain top shape matches. 


TABLE IV. Retrieval rate(Bulls eye Score) oe MPEG-7 

DATASET- A COMPARATIVE ANALYSIS 


DataSet 

MPEG-7 

css m — 

75.44 

Visual parts (6| 

76.45 

sc H 

76.51 

CPDH "p2| 

76.56 

Aligning curves jisj 

78.16 

Proposed Approach 

79.38 

ssc 

79.92 


VII. Experimental Results and Discussions 

In this section, we present the experimental results con¬ 
ducted on the standard shape datasets namely: Kimia-99, 
Kimia-216 and MPEG-7. The performance of the proposed 
approach is demonstrated through Bull’s eye score |[^. All 
experiments are conducted using MATLAB tool and tested on 
Pentium(R) dual core CPU with 3GB RAM on Windows-7. 


have also given the retrieval results obtained due to Belongie 
et al. j^and Shu and Wu 1^ . One can notice here too that the 
proposed approach works on par with the existing methods. 


□ n II s 


A. Experimental results on Kimia-99 

Kimia-99 shape dataset consists of 9 classes with 11 
samples in each class. The 99 shapes belong to 9 different 
classes are shown in Figure The retrieval rate depicting 
the top 10 matching shapes is presented in Table |I| We have 
also given the retrieval results obtained due to Belongie et 
al. m, Shu and Wu ^ and Tu and Yuille ||2^. It shall 
be observed from Table]^ that the proposed approach exhibits 
better performance when compared to other approaches. 
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Fig. 5. Kimia 99 dataset: 9 classes with 11 samples in each class 


B. Experimental results on Kimia-216 

Kimia-216 shape dataset consists of 18 classes with 12 
samples in each class. One such sample from each class is 
shown in Figure The top 11 closest matches obtained due to 
the proposed methodology is shown in Table In Table we 


Fig. 6. Kimia 216 dataset : A sample shape from each class 


C. Experimental results on MPEG-7 : 


MPEG-7 dataset consists of 1400 shapes from 70 classes 
with each class consisting of 20 shape samples. Figureshows 
one sample shape from each class of MPEG-7. The dataset is 
challenging due to the presence of shapes that are visually 
dissimilar from other members of their class and shapes that 
are highly similar to members of other classes. The retrieval 
rate of the proposed approach on MPEG-7 is calculated and 
the results obtained due to proposed approach along with the 
results of the state-of-art algorithms is tabulated in Table EYl 
The retrieval rate obtained due to proposed approach for each 
shape class considering top 40 retrievals is shown in Figure 
and also tabulated in Table |II^ The retrieval rate depicting top 
12 closest matching shape is given Table [Vl| 


Table [V] provides overall performance of the proposed 
approach measured using bull’s eye score on Kimia-99, Kimia- 
216 and MPEG-7 shape datasets. 
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Fig. 7. MPEG-7 dataset: A sample shape from each class 
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Fig. 8. Clas-wise retrieval results for MPEG-7 dataset 


TABLE III. Clas-wise retrieval results eor MPEG-7 dataset 


Class 

Device9 

Device6 

Bird 

Horse 

Octopus 

Chicken 

Classic 

Ray 

Bell 

device 

MPEG-7 

5 

6 

7 

7 

8 

9 

9 

9 

10 

10 

Class 

Device2 

Turtle 

Cattle 

Device 1 

elephant 

Imfish 

Tree 

Device3 

Device4 

Prog 

MPEG-7 

10 

10 

11 

11 

11 

11 

11 

12 

12 

12 

Class 

SeaSnake 

Watch 

Devices 

Pencil 

deer 

spoon 

Spring 

Device7 

Jar 

Pocket, 

MPEG-7 

12 

12 

13 

13 

14 

14 

14 

15 

15 

15 

Class 

Butterfly 

Camel 

Lizzard 

Beetle 

Device8 

Ply 

Pork 

Hat 

Horseshoe 

Bat 

MPEG-7 

16 

16 

16 

17 

17 

17 

17 

17 

17 

18 

Class 

Dog 

Hammer 

Misk 

Crown 

Pish 

Guitar 

PersonalCar 

bone 

Comma 

Glas 

MPEG-7 

18 

18 

19 

19 

19 

19 

19 

20 

20 

20 

Class 

HCircle 

Heart 

Apple 

Bottle 

Chopper 

Cup 

Brick 

Car 

Carriage 

CellularPhone 

MPEG-7 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

Class 

Children 

Chopper 

Elatfish 

Eountain 

Key 

Rat 

Shoe 

Stef 

Teddy 

Truck 

MPEG-7 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 


TABLE V. Bull Eye Score(Retrieval Rate) oe proposed 
METHOD EOR KiMIA 99, KiMIA 216, MPEG-7 DATASET 


DataSet 

Bull eye test 

Kimia-99 

99.08 

Kimia-216 

95.25 

MPEG-7 

79.38 


TABLE VI. Top 12 closest matching shapes 


Dataset 

1st 

2nd 

3rd 

4th 

5th 

6th 

7th 

8th 

9th 

lOth 

llfl] 

12tl 

MPEG 

-7 

1400 

1345 

1276 

1221 

1157 

1113 

107( 

1022 

995 

961 

933 

898 


VIIL Conclusion 

In this paper, we explored the application of morpho¬ 
logical pattern spectra and Local binary pattern method on 
binary shapes to extract the dominant features based on 
pattern spectra followed by classification using EMD. The 


Pattern spectra is based on the concept introduced by Maragos 
transformation(Thickness map). The relationship between the 
Maragos morphological pattern spectrum and the histogram 
is well utilized here. Also the fusion of LBP method along 
with the pattern spectra gives better results at a comparable 
time complexity. It is proved that the morphological spectra 
and LBP are invariant to rotation and shift of the shape. 
This method is proven to be computationally efficient. It is 
experimentally shown that morphological spectra together with 
LBP, in combination with EMD metric, with LI as the ground 
distance can be successfully used in the shape based object 
classification. The proposed approach is relatively simple to 
implement, and accurate in terms of shape representation and 
classification. 
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